试除法判定约数(sqrt(n))
void get_divisors(int n)
{
    vector<int> res;
    for (int i = 1; i <= n / i; i++) {
        if (n % i == 0) {
            res.push_back(i);
            if (i != n / i) {  // 避免 i==n/i, 重复放入 （n是完全平方数
                res.push_back(n / i);
            }
        }
    }
}

基本思想：
如果 N=p1c1*p2c2*…*pkck
约数个数：(c1+1)*(c2+1)*..(ck+1)
约数之和：f(n)=(p1^0+p1^1+p1^2+…p1^a1)(p2^0+p2^1+p2^2+…p2^a2)…(pk^0+pk^1+pk^2+…pk^ak）

最大公约数,辗转相减优化到辗转相除
int gcd(int a,int b)
{
    return b==0 ? a:gcd(b,a%b);
}
最小公倍数 lcm(a,b) * gcd(a,b) = a * b

由算数基本定理，将两个数分解为质数乘积形式，
其最大公约数为取每个质数的较小次幂的乘积，最小公倍数为取每个质数的较大次幂的乘积